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We estimate the neutrino emission from individual gamma-ray bursts ob- 
served by the BATSE detector on the Compton Gamma-Ray Observatory. Neu- 
trinos are produced by photoproduction of pions when protons interact with pho- 
tons in the region where the kinetic energy of the relativistic fireball is dissipated 
allowing the acceleration of electrons and protons. We also consider models where 
neutrinos are predominantly produced on the radiation surrounding the newly 
formed black hole. From the observed redshift and photon flux of each individual 
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burst, we compute the neutrino flux in a variety of models based on the assump- 
. tion that equal kinetic energy is dissipated into electrons and protons. Where not 

measured, the redshift is estimated by other methods. Unlike previous calcula- 
tions of the universal diffuse neutrino flux produced by all gamma-ray bursts, the 
individual fluxes (compiled at http://www.arcetri.astr0.1t/~ dafne/grb/) can be 
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directly compared with coincident observations by the AMANDA telescope at 
the South Pole. Because of its large statistics, our predictions are likely to be 



representative for future observations with larger neutrino telescopes. 



Subject headings: gamma rays: bursts — acceleration of particles — neutrinos 



1. Introduction 

The leading models for Gamma- Ray Bursts (GRB), bursts of 0.1-1 MeV photons typ- 
ically lasting for 0.1-100 seconds (Fishman & Meegan 1995), involve a relativistic wind 
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emanating from a compact central source. The ultimate energy source is rapid accretion 
onto a newly formed stellar mass black hole. 

Observations suggest that the prompt 7-ray emission is produced by the dissipation 
(perhaps due to internal shocks) of the kinetic energy of a relativistically expanding wind, 
i.e. a "fireball" . Both synchrotron and inverse Compton emissions from the shock-accelerated 
electrons have been proposed as the GRB emission mechanism. 

In this paper, we study in detail the production of neutrinos by protons accelerated 
along with electrons. We assume that equal energy of the fireball is dissipated in protons 
and electrons (or photons). This is the case in models where GRBs are the sources of the 
highest energy cosmic rays. The basic idea is that the protons produce pions decaying into 
neutrinos in interactions with the fireball photons, or with external photons surrounding the 
newly formed black hole. 

Where previous calculations have estimated the universal diffuse flux of neutrinos pro- 
duced by all GRBs over cosmological time, we estimate the flux from individual GRBs ob- 
served by the BATSE (Burst And Transient Source Experiment) experiment on the Compton 
Gamma-Ray observatory. The prediction can be directly compared with coincident obser- 
vations performed with the AMANDA detector. Having these observations in mind, we 
specialize on neutrino emission coincident in time with the GRB. Opportunities for neutrino 
production exist after and, in some models, before the burst of 7-rays, e.g. when the fireball 
expands through the opaque ejecta of a supernova. 

The calculations are performed in two models chosen to be representative and rather 
different versions of a large range of competing models. The first is generic for models where 
an initial event, such as a merger of compact objects or the instant collapse of a massive star 
to a black hole, produces the fireball (Waxman & Bahcall 1997 ;Guetta, Spada & Waxman 
2001a, GSW hereafter). We calculate the neutrino production by photomeson interactions of 
relativistic protons accelerated in the internal shocks and the synchrotron photons that are 
emitted in these shocks. The neutrinos produced via this mechanism have typical energies 
of ~ 10 14 — 10 15 eV, and are emitted in coincidence with the GRBs, with their spectrum 
tracing the GRB photon spectrum. For an alternative model, we have chosen the supranova 
model where a massive star collapses to a neutron star with mass ~ 2.5 — 3M , which 
loses its rotational energy on a time scale t s d of weeks to years, before collapsing to a black 
hole, thus triggering the GRB. Following Guetta and Granot (Guetta & Granot 2002a), we 
calculate the neutrino flux from interactions of the fireball protons with external photons in 
the rich radiation field created during the spindown of the supra-massive pulsar. Production 
on external photons turns out to be dominant for a wide range of parameters, t s d < 0.2 yr for 
a typical GRB and t s< \ < 2 yr for X-ray flashes. The neutrinos produced via this mechanism 
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have energies e v ~ 10 15 — 10 17 (10 19 ) eV for typical GRB (X-ray flashes) and are emitted 
simultaneously with the prompt 7-ray (X-ray) emission. Their energy spectrum consists of 
several power law segments and its overall shape depends on the model parameters, especially 
t s d- If the mass of the supernova remnant is of the order of ~ 1OM , and if the supernova 
remnant shell is clumpy, then for time separation t s d < 0.1 yr the SNR shell has a Thomson 
optical depth larger than unity and obscures the radiation emitted by the GRB. Therefore, 
for t s d ^0.1 yr, the z/'s would not be accompanied by a detectable GRB providing us with 
an example of neutrino emission not coinciding with a GRB display. 

As previously mentioned, to realistically estimate the neutrino fluxes associated with 
either of these models, we turn to the GRB data collected by BATSE. The BATSE records 
include spectral and temporal information which can be used to estimate neutrino spectra 
for individual bursts. We will perform these calculations for two version of each of the 
two models previously described. The two versions correspond to alternative choices of 
important parameters. The wide variety of GRB spectra results, not surprisingly, in a wide 
range of neutrino spectra and event rates. For approximately 800 bursts in the BATSE 
catalog, and for four choices of models, we have calculated the neutrino spectra and the 
event rates, coincident with GRBs, for a generic neutrino telescope. With 800 bursts, the 
sample should also be representative for data expected from much larger next-generation 
neutrino observatories. 

Neutrino telescopes can leverage the directional and time information provided by 
BATSE to do an essentially background-free search for neutrinos from GRBs. Individual 
neutrino events within the BATSE time and angular window are a meaningful observation. 
A generic detector with 1km 2 effective telescope area, during one year, should be able to 
observe 1000 bursts over Arc steradians. Using the BATSE GRBs as a template, we predict 
order 10 events, muons or showers, for both models. The rates in the supranova model 
depend strongly on t sd . In this model we anticipate ~ 7 events per year assuming t sd — 0.07 
yr, but only one event per ten years for t s d — 0.4 yr. We will present detailed tabulated 
predictions further on. They can be accessed at http://www.arcetri.astro.it/~dafne/grb/. 

Short duration GRBs, characterized by lower average fluences, are less likely to produce 
observable neutrino fluxes. 

We find that GRBs with lower peak energies, X-ray flash candidates, yield the largest 
rates in the supranova model. For instance, for t sd ~ 0.07 yr, we predict one event (muon 
or shower) for every 1000 bursts. If only 100, or so, X-ray flashes occur per year, as obser- 
vations suggest, this will be difficult to observe. However, such events may be considerably 
more common and may contribute significantly to the diffuse high energy neutrino flux. Ob- 
servations of neutrinos from this class of GRBs would be strong evidence for a supranova 
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progenitor model. 

The AMANDA collaboration has collected neutrino data in coincidence with BATSE 
observations (Barouch & Hardtke 2001). It operated the detector for 3 years (1997-1999) 
with an effective area of approximately 5,000 m 2 for GRBs events. In early 2000 the expanded 
detector reached an effective area of roughly 50,000 m 2 . Unfortunately, only ~ 100 coincident 
bursts could be observed with the completed detector before BATSE operations ceased in 
June 2000. With effective areas significantly below the canonical square kilometer discussed 
in this paper, AMANDA is not large enough to test the GRB models considered here. We 
estimate only 0.08 events in 1997-1999 and 0.3 in 2000. 

Our estimate of 1-10 events in 1 year for a telescope with 1 kilometer square telescope 
area is consistent with previous, burst-averaged, determinations of the GRB neutrino flux 
(Alvarez, Halzen & Hooper 2000, Dermer & Atoyan 2003). Note that the effective area 
of IceCube for GRB will significantly exceed this reference value (PDD). The rate can be 
understood by a back-of-the-envelope estimate. A typical GRB produces a photon fluence 
on the order of 10~ 5 ergs/cm 2 , which we assume to be equal to the energy in protons. If 
20% of the proton energy is converted into pions, half to charged pions, and one quarter 
of the charged pion energy to muon neutrinos in the n — > /i — > e decay chain, then ~ 
5 x 10~ 7 ergs/ cm 2 /E u neutrinos are generated. For a typical neutrino energy of E v ~ 100 TeV, 
this yields ~ 30 neutrinos per square kilometer. At this energy the probability that a 
neutrino converts to a muon within range of the detector is 10~ 4 (Gaisser, Halzen, Stanev 
1995). Therefore 3 x 10~ 3 muons are detected in association with a single GRB. With over 
1000 GRBs in one year, we estimate a few muons per year in a kilometer-square detector. 
Fluctuations in fluence and other burst characteristics enhance this estimate significantly 
(see Figure 2, for instance), however, absorption of neutrinos in the Earth can reduce this 
number. These effects are included in the calculations of actual event rates throughout this 
paper. 

Our estimates, while observable with future kilometer-scale observatories, may be con- 
servative. We already mentioned bursts with no counterpart in photons. Also, occasional 
nearby bursts, much like supernova, could exceptionally provide higher event rates than our 
calculations reflect. We would like to point out the fact that the aim of the paper is not 
to do more precise calculations than the ones already done in the literature. We just want 
to generate results that can be compared with experiments that do coincidents observations 
with satellites. 

The outline of the paper is somewhat unconventional. The detailed results are collected 
in a data archive at http:/ /www. arcetri.astro.it/~dafne/grb. The details of the calculations 
of the neutrino fluxes in the two models are described in appendices A and B. In appendix 
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C we collected the methods used to evaluate the rates of neutrino-induced muons, taus and 
showers in a generic detector with a given effective telescope area. The main body of the 
paper is organized as follows. In §2, we review GRB models and describe the mechanisms 
for neutrino production in GRBs. In §3 and §4 we describe the BATSE catalog of GRBs and 
we identify subclasses of events: long duration GRBs with and without measured redshift, 
short duration GRBs and X-ray flash candidates. In §5 we summarize our simulation of the 
response of a generic high energy neutrino telescope to the predicted neutrino fluxes. In §6 
we analyse some anomalous bursts. Results and conclusions are collected in §7. 



2. GRB Models and Neutrino Emission 

Progenitor models of GRBs are divided into two main categories. The first category 
involves the merger of a binary system of compact objects, such as a double neutron star 
(Eichler et al. 1989), a neutron star and a black hole (Narayan, Pacyhski & Piran 1992) or 
a black hole and a Helium star or a white dwarf (Fryer & Woosley 1998; Fryer, Woosley & 
Hartmann 1999). The second category involves the death of a massive star. It includes the 
failed supernova (Woosley 1993) or hypernova (Pacyhski 1998) models, where a black hole 
is created promptly, and a large accretion rate from a surrounding accretion disk (or torus) 
feeds a strong relativistic jet in the polar regions. This type of model is known as the collapsar 
model. An alternative model within this second category is the supranova model (Vietri & 
Stella 1998), where a massive star explodes in a supernova and leaves behind a supra- massive 
neutron star (SMNS), of mass ~ 2.5 — 3M . It subsequently loses its rotational energy on 
a time scale t sd of order weeks to years until it collapses to a black hole. This triggers the 
GRB. Long GRBs (with a duration > 2 s) are usually attributed to the second category of 
progenitors, while short GRBs are attributed to the first category. 

We select two models to investigate the opportunities for neutrino production in a GRB. 
The analysis can be easily extended to other models (see for example Razzaque, Meszaros 
& Waxman 2002; Dermer & Atoyan 2003). The first model is based on the standard fireball 
phenomenology where electrons and protons are shock accelerated in the fireball. Pions and 
neutrinos are produced by photoproduction interactions when the protons coexist in the 
fireball with photons. These are produced by synchrotron radiation of accelerated electrons. 
For a second model we turn to the supranova scenario where the supra-massive pulsar loses 
its rotational energy through a strong pulsar wind. This pulsar wind creates a rich external 
radiation field before the collapse to the final black hole and the creation of the GRB fireball. 
It is referred to as the pulsar wind bubble (PWB) (Konigl & Granot 2002; Inoue, Guetta 
& Pacini 2002; Guetta & Granot 2003) and provides a target for the photoproduction of 
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neutrinos by fireball protons. Let's note in passing that the supranova model has several 
advantages compared to other collapsar models: (i) the jet does not have to penetrate the 
stellar envelope (Vietri & Stella 1998), (ii) it can naturally explain the X-ray line features 
observed in several afterglows (Piro et al. 2000; Lazzati, et al. 2001; Vietri et al. 2001), and 
(hi) the large fraction of the internal energy in the magnetic field and in electrons observed 
in the afterglow emission arise naturally (Konigl & Granot 2002; Guetta & Granot 2003). 

In all of the different scenarios mentioned above, the final stage of the process consists of 
a newly formed black hole with a large accretion rate from a surrounding torus, and involve 
a similar energy budget (< 10 54 ergs). Observations suggest that prompt 7-ray emission is 
produced by the dissipation of the kinetic energy within the fireball, due to internal shocks 
within the flow that arise from variability of the Lorentz factor, T, on a time scale t v . The 
afterglow emission is thought to arise from an external shock that is driven into the ambient 
medium as it decelerates the ejected matter (Rees & Meszaros 1994; Sari & Piran 1997). In 
this so called 'internal-external' shock model, the duration of the prompt GRB is directly 
related to the time during which the central source is active. The emission mechanism is 
successfully described by synchrotron radiation from relativistic electrons that radiate in the 
strong magnetic fields. These are close to equipartition values within the shocked plasma. 
An additional radiation mechanism that may also play some role is synchrotron self-Compton 
(SSC) (Guetta & Granot 2002b), which is the upscattering of the synchrotron photons by 
relativistic electrons to higher energy. 

Protons are expected to be accelerated along with the electrons in the region where 
the wind kinetic energy is converted into internal energy due to a dissipation mechanism 
like internal shocks. The conditions in the dissipation region allow proton acceleration up 
to £ p ,max ^ 10 20 eV (Waxman 1995; Vietri 1995). The energy in 7-rays reflect the fireball 
energy in accelerated electrons and afterglow observations indicate that accelerated electrons 
and protons carry similar energy (Freedman & Waxman 2000). Our basic assumption in 
calculating neutrino emission from GRBs is that equal amounts of energy go into protons 
and photons. In models where GRB protons are the source of the highest energy cosmic 
rays, this assumption is supported by the approximate equality of the 7-ray fluence of all 
GRBs and the total energy in extragalactic cosmic rays. 

Both internal shocks, responsible for the prompt GRB emission, and the external shock, 
responsible for the afterglow emission, have been proposed as possible sources of the highest 
energy cosmic rays (Waxman 1995 and Vietri 1995, respectively). A comparison between the 
two mechanisms has been done by Vietri, De Marco and Guetta (2003), but it is not easy 
to conceive, at this point, an observational test capable of distinguishing between the two 
models. The only hope appears to observe the production of high energy neutrinos which 
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must accompany the in situ acceleration of particles. Occasionally, ultra-high energy protons 
will produce pions and neutrinos in collisions with photons in photon rich environment 
provided by the post-shock shells or by the pulsar wind bubble. If the protons are accelerated 
in internal shocks, the neutrinos produced will arrive at Earth simultaneously with the 
photons of the burst proper and will have an energy ~ 10 15 — 10 16 eV (Waxman and Bachall 
1997; GSW; Guetta & Granot 2002a). If accelerated in external shocks, they will arrive at 
Earth simultaneously with the photons of the afterglow and will have a higher energy, > 10 17 
eV (Vietri 1998a, 1998b, Waxman & Bahcall 2000). 

General phenomenological considerations indicate that gamma-ray bursts are produced 
by the dissipation of the kinetic energy of a relativistic expanding fireball. Internal shocks 
that are mildly relativistic are believed to dissipate the energy. Therefore, the proton energy 
distribution should be close to that for Fermi acceleration in a Newtonian (non-relativistic) 
shock, dn p /de p oc e p 2 . Moreover, the power law index of the electron and proton energy 
distributions are expected to be the same, and the values inferred for the electron distribution 
from the observed photon spectrum are dn e /de e oc e~ p with p ~ 2 — 2.5. We shall, therefore, 
adopt drip/dEp oc e p 2 . Plasma parameters in the dissipation region allow proton acceleration 
to energies greater than 10 20 eV (Waxman 1995, Vietri 1995). 

We will assume that the fireball is spherically symmetric. Note, however, that a jet-like 
fireball behaves as if it were a conical section of a spherical fireball as long as 9j > 1 /T, where 
9j is the jet opening angle and T > 300 is the wind Lorentz factor. Therefore, our results 
apply without modification to a jet-like fireball. For a jet-like wind, the luminosity, L, in 
our equations should be understood as the luminosity of the fireball inferred by assuming 
spherical symmetry. 

We have relegated all details of the calculation of neutrino production via photomeson 
interaction with GRB photons and PWB photons to appendices A and B, respectively. 
Throughout the paper, we will refer to the models by the following convention: 

(i) Model 1: Neutrino flux from the interaction of high-energy protons with GRB photons. 
The fraction of proton energy transfered to pion energy is set to 0.2 as indicated by 
simulations of GSW (see apppendix A). 

(ii) Model 2: Neutrino flux from the interaction of high-energy protons with GRB photons. 
The fraction of proton energy transfered to pion energy is calculated as described in 
appendix A (see appendix A). 

(iii) Model 3: Neutrino flux from the interaction of high-energy protons with PWB photons, 
as in the supranova progenitor model. The time scale between supernova and GRB is 
set to t sd = 0.4 yr (see appendix B). 
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(iv) Model 4: Neutrino flux from the interaction of high-energy protons with PWB photons, 
as in the supranova progenitor model. The time scale between supernova and GRB is 
set to t sd = 0.07 yr (see appendix B). 

3. The BATSE Catalog 

BATSE, the Burst And Transient Source Experiment, was a high energy astrophysics 
experiment launched on the Compton Gamma- Ray Observatory in 1991. BATSE, between 
its launch and the termination of its orbit in 2000, has observed and recorded data from over 
8000 events including gamma-ray bursts, pulsars, terrestrial gamma-ray flashes, soft gamma 
repeaters and black holes. 

The data BATSE recorded from gamma-ray bursts is publically available in the cur- 
rent BATSE catalog at http://f64.nsstc.nasa.gov/batse/grb/catalog/. For a description see 
Paciesas et al. 1999. The catalog includes information on the spectrum, time and location 
of each triggered burst. Each triggered event has been assigned a BATSE trigger number 
(between 105 and 8121 for GRBs) which we use to identify individual bursts. 

Spectral information is recorded in four energy channels, 20-50 keV, 50-100 keV, 100-300 
keV and above 300 keV. Using these four fluence measurements, we have fitted the spectrum 
of each GRB to a broken power law, treating the break energy, both spectral slopes and 
the normalization as free parameters; see Eq.(A6). We determine the Lorentz factor of the 
relativistic expanding ejecta using the break energy through Eq.(Al8). For bursts with an 
observed break energy above 300 keV, or below 50 keV, it is difficult to determine both the 
break energy and the power law slope. For high energy breaks, the impact on this ambiguity 
is not critical. As explained in appendix A, the Lorentz factor is not dependent on the fit 
because it is fixed by the observed high energy of the event and the requirement that the 
fireball be optically thin. For very low spectral breaks (below 50 keV), for instance in events 
we classify as X-ray flash candidates, we acknowledge a significant degree of uncertainty in 
the Lorentz factor calculation and resulting neutrino spectra. In this case, the spectral break 
is only uncertain to about a factor of 2 or 3. 

Detailed temporal information is available in the BATSE catalog, in the form of light 
curves. The BATSE time resolution varies between 2.048 seconds and 0.016 seconds. A 
resolution of 0.064 seconds is available for all bursts during the time following the trigger. 

In the framework of the internal shock model, we need a variability time, t v , < 0.01s in 
order to get the lMeV 7 — rays (Guetta, Spada & Waxman 2001b, Waxman 2001). In fact 
with a value of the Lorentz factor larger than the minimal value needed to be optically thin 
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up to 100 Mev (~ 250 see Eq.Al7) the variability time, t v , has to be < 10ms to get the 1 
MeV 7-rays. 

How well the data support the model is still controversial and a detailed analysis on 
this issue is out of the aim of our paper. Since we refer to this model for our analysis we 
consider a value of t v = 10 ms for long duration GRBs in the rest of the paper. For short 
duration bursts we take t v — 0.1 ms and 0.050 seconds for X-ray flash candidates. 



We have divided the list of BATSE GRBs into four different classes: 1) long duration 
bursts (duration > 2 s) with measured redshfit, 2) long duration bursts without measured 
redshift, 3) short duration bursts and 4) X-ray flashes. Major differences in temporal and 
spectral properties of long duration GRBs, short duration GRBs and X-ray flashes has lead 
to some speculation that they may involve different progenitors or mechanisms. 



By observing the optical afterglow of GRBs, it is possible to measure spectral lines and, 
therefore, the redshift of an individual burst. Although to date the X-ray afterglow of on the 
order of 100 long duration GRBs have been observed, only 31 of them have been observed in 
the optical making a determination of the redshift possible. (No redshift has been identified 
for short GRBs). We first consider 13 of these that have a complete BATSE record. 

The relationship between the comoving distance to an object and its redshift is given 



Where Q\ — 0.7 and Q M ~ 0.3 are the fractions of the critical density of the Universe in 
dark energy and matter, respectively. H is Hubble's constant. 

Once the distance to a GRB is known, and its fluence (or flux) has been measured, 
the isotropic-equivalent 7-ray luminosity can be calculated. From this and the value of 
t v = 0.01 s, together with the break energy of the GRB photon spectrum, we estimate 
the bulk Lorentz factor using Eqs.(Al7, A18). Setting the efficiency for pion production 
f w = 0.2 for model 1 and using Eq.(Al6) for model 2, we estimate the fraction of proton 
energy transfered into pions. Using Eqs.(Al0, All) the energy scale of synchrotron losses, 



4. GRB Classes 



4.1. Long Duration GRBs With Measured Redshift 



by: 




(1) 
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sl, is determined. From these informations we can determine the neutrino flux at Earth for 
each of the four models. For a more detailed discussion see appendices A and B. 



4.2. Long Duration GRBs Without Measured Redshift 

For the majority of the long duration bursts in the BATSE catalog no redshift is avail- 
able. In this case, the 7-ray luminosity of long duration bursts, and, therefore, the distance, 
can be estimated by assuming a relationship between the observed variability and the lumi- 
nosity of a GRB (Lloyd-Ronning & Ramierz-Ruiz 2002; Zhang & Meszaros 2002; Kobayashi, 
Ryde & MacFadyen 2002). Note that the variability of a GRB is not the same as its vari- 
ability time, t v , previously introduced. 

The variability of burst is a measure of fluctuations in the temporal structure of the 
burst. It is defined such that pure noise should have a variability of zero, while the most 
variable bursts have very sudden and distinctive temporal features. We use the following 
definition of variability (Fenimore & Ramirez- Ruiz 2000): 



va-0.24 1 ^ ( F i - -Fo.30xr 90 ) - (Fj + B) 



where Y = (l + z)/3, Fi is the background subtracted fluence in a time bin i, B is the average 
background in a single time bin, F peak is the maximum fluence and F 30xTgo is the average 
fluence over a time period centered at time bin % of length 30% of the T 90 time (duration) of 
the burst. 

The GRBs with observed redshifts have been used to empirically derive a relationship 
between the variability and the luminosity of the burst (Fenimore & Ramirez- Ruiz 2000; 
Reichart et al. 2001, Reichart & Lamb 2001): 

L/dVL = 3.1 x 10 5f V 3 - 35 erg s _1 . (3) 

A relationship between luminosity and the time lag between the peaks for light curves in 
different energy bands, has been observed in the GRB redshift data (Norris, Marani & 
Bonnell 2000), but appears to be less reliable. Therefore, we will only consider the luminosity- 
variability relationship. 

Together with the fluence (flux) of a burst, the distance of a burst can be determined 
by the luminosity. Note that the variability, and therefore luminosity, of a burst depends on 
the redshift or distance to the burst. Therefore, we must do this calculation by iteration. 
The end result is a value of the luminosity and redshift for each burst. We would like to 
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enphasize the fact that there are a lot of uncertainties in this way to estimate the redshift, 
however the knoweledge of the redshift is not so important in our analysis since the minimum 
variability time scales and the fluence are the most important quantities. 

It is difficult to reliably calculate the variability of bursts with low flux. For this reason, 
we only consider bursts with a peak flux greater than 1.5 photons/cm 2 sec (over a 0.256 s 
time scale) and with at least 30 time bins (of 0.064 s width) of 5 sigma or more above the 
average background. The GRBs which do not meet these requirements have low fluence and 
are therefore likely to yield a low neutrino flux anyway. After these criteria were applied, 
566 long duration bursts without measured redshift are left, making it our largest class. 

Even for bursts which meet the above criteria, the luminosity estimated is only accurate 
to an order of magnitude. This corresponds to uncertainties of a factor of 2 or 3 in the 
fraction of proton energy transfered into pions and in the synchrotron loss energy. 

4.3. Short Duration GRBs 

Short duration GRBs, with no observation of an optical afterglow and, therefore, no 
measurement of redshift, cannot have a relationship between variability and luminosity em- 
pirically established. Additionally, variability is difficult to measure for short duration bursts. 
Left with no way to measure the 7-ray luminosity of, or distance to, a short duration GRBs, 
we choose to set z — 1 for each burst. This introduces greater uncertainty than in the long 
duration GRBs case but, given the ambiguities in the burst characteristics, it is the best that 
can be done at this time. To be consistent, we considered only short bursts with a peak flux 
greater than 1.5 photons/cm 2 sec (over a 0.256 s time scale), as we did with for long duration 
GRBs. After this criteria was applied, 199 short duration bursts remained in this category. 

Temporal structure and variations appear to occur on shorter time scales for short 
compared to long GRBs (McBreen et al. 2002). We therefore use a time scale of temporal 
fluctuations of t v = 0.001 s for all short duration bursts, as opposed to the value of t v = 0.01 s 
used for long duration bursts. It is also interesting to note that short duration bursts 
generally have a somewhat harder spectrum and higher peak energy than long duration 
GRBs (Paciesas et al. 2001). 

4.4. X-Ray Flash Candidates 

X-ray flashes are a newly discovered class of fast transient sources. The BeppoSAX 
experiment's wide field cameras have observed such events at a rate of about four per year 
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(Heise et al. 2001, Kippen et al. 2002), implying a total rate on the order of 100 per 
year. These events typically have peak energies as low as 2-10 keV and durations of 10-100 
seconds. The BeppoSAX experiment discriminates X-ray flashes from standard GRBs by the 
non-detection of a signal above 40 keV with the BeppoSAX GRB monitor. More generally, a 
large ratio of X-ray to 7-ray fruence is the differentiating characteristic of X-ray flashes from 
GRBs. It has been suggested, however, that X-ray flashes are a low peak energy extension 
of gamma-ray bursts (Heise et al. 2001, Kippen et al. 2002). 

The final class considered here consists of BATSE events which may be X-ray flashes. 
We identify 15 events in this class with spectra that peak below 50 keV, although it is 
difficult to determine accurately where the peak occurs because the sensitivity of BATSE is 
somewhat poor in this energy range. These are long duration bursts, and typically have a 
hard spectrum; several have (3 in Eq.(A6) larger than 1.5, and no observed flux in BATSE's 
fourth energy channel (above 300 keV). 

Again, with no measured redshift, and limited temporal information, we cannot deduce 
the luminosities of these events. We set z — 1 for each event and calculate its luminosity 
accordingly. It is thought that the radius of collisions in X-ray flashes is typically larger than 
in other GRBs and, therefore, the time scale of fluctuations will generally be larger (Guetta, 
Spada & Waxman 2001b). We therefore choose t v = 0.05 s. For some X-ray flashes the peak 
energy is very low (20-30 keV) and the Lorentz factor accordingly very high. Increasing 
the time scale of fluctuations has the additional effect of lowering the Lorentz factor to a 
reasonable value in these extreme cases (see Eq.(Al8)). 

For X-ray flashes, with a very large Lorentz factor and a long time scale of fluctuations, 
we expect that a very small fraction of proton energy converted to pions; see Eqs.(A15,A16). 
We therefore expect low neutrino fluxes from proton interactions with GRB photons. For 
the models involving proton interactions with photons in a surrounding pulsar wind bubble, 
however, the rates can be quite high (Guetta & Ganot 2002a). 

There is another class of objects that are the non-triggered bursts. They are found in 
the BATSE data in off-line analysis; see e.g. Stern and Tikhomirova 

(http://www.astro.su.se/groups/head/grb_archive.html). They are not energetic enough to 
trigger in real time. These can be GRBs with low kinetic luminosity and will be very weak 
neutrino sources for all the four models, therefore we have decided to neglect them in our 
analysis. 
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4.5. Neutrino GRB Database 

We have constructed a database publicly accessible online with a complete list of all 
of the GRB characteristics and associated neutrino event rates for the approximately 800 
bursts we have considered in this analysis. The data was populated in a MySQL database, 
interacting through a user- friendly web interface developed with PERL. The database is 
searchable by GRB number and/or by GRB class, and is capable of searching for only 
bursts included in the AMANDA analysis or all bursts in this work. The rates are for a 
generic neutrino telescope provided the threshold is sufficiently low for observing the neutrino 
fluxes predicted. How we transform the GRB neutrino fluxes into observed event rates is 
the topic of the next section. Representative results are tabulated in tables 1 through 12. 
The database also contains the predicted neutrino spectra for each burst; these can be 
directly combined with the simulation of a specific detector. The database is accessible at 
http: / / www.arcetri.astro.it / ~dafne / grb. 

5. Neutrino Telescopes and Event Simulation 

Large volume neutrino telescopes are required to observe and measure the neutrino flux 
from GRBs. Current experiments, such as AMANDA (Andres et al. 2001) at the South 
Pole, or future (Aslanides et al. 1999) and next generation experiments, including IceCube 
(see http://icecube.wisc.edu/) with a full cubic kilometer of instrumented detector volume, 
i.e. with > 1km 2 telescope area, are designed to observe high energy cosmic neutrinos with 
energies expected from GRBs. 

Neutrino telescopes detect the Cherenkov light radiated by showers (hadronic and elec- 
tromagnetic), muons, and taus that are produced in the interactions of neutrinos inside or 
near the detector. Muons are of particular interest because, at the energies typical for GRB 
neutrinos, they travel kilometers before losing energy. The dominant signal is, therefore, 
through-going muons. IceCube can measure the energy and direction of any observed muon. 
The angular resolution is less than 1° — 3° while the energy resolution is approximately a 
factor of three. Signal and background muons may, therefore, be differentiated with a simple 
energy cut. For shower events the energy measurement improves significantly, being better 
than 20%, but reconstructing their direction is challenging, the angular resolution being of 
order 10 degrees. 

For high energies, when tau decay is sufficiently time dilated, taus have a range similar 
to or larger than muons, and so the dominant tau signal is from through-going taus. These 
events have a characteristic signature consisting of a "clean" minimum-ionizing track despite 
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its long range inside the detector. We will therefore consider all events in which a tau 
track passes through the detector in the direction and at the time of a GRB. We assume 
that taus and muons are distinguishable at all energies when specializing to rare events in 
coincidence in time and direction with a GRB. We realize that, in general, it may be difficult 
to distinguish a muon of energy < 200 GeV, that is expected to lose relatively little energy 
from catastrophic processes, from a very high energy tau. Tau signatures however become 
dramatic when they decay inside the detector (lollipop events), or when the tau neutrino 
interacts and the produced tau decays into showers inside the detector (double bang events). 

To evaluate the prospects for GRB neutrino observations, it is essential to determine 
the rate of muon, shower and tau events. These calculations are each described in appendix 
C. For a review of high energy neutrino astronomy; see (Halzen & Hooper 2002; Learned & 
Mannheim 2000). 

6. Anomalous Bursts 

There are a few GRBs we have considered in this analysis which are anomalous for a 
variety of reasons. We briefly mention these in this section. 

(i) GRB 6707 is a burst with a measured redshift of 0.0085, yet a relatively low ffuence 
of 4 x 10~ 6 erg/cm 2 . Together, this implies a luminosity on the order of 10 45 erg/s, 
well below the normal range considered. Our calculation of the Lorentz factor, which 
depends on the luminosity of the GRB, yields a value of about 25,000, much larger 
than the normally allowed range. Results in model 2 are, therefore, not particularly 
realistic for this particular burst (the rates are actually very low for this model). The 
other models are only affected by this in the calculation of the synchrotron energy loss 
scale. 

(ii) GRBs 7648, 6891 and 1997 each have Lorentz factors below 100. In all three cases, 
the burst has been found to have a low luminosity, ~ 10 49 erg/s, which contributes 
to this result. The results of models 2, 3 or 4 are only affected by this in their low 
synchrotron energy loss scales, and are, therefore, conservative. Results for model 1 
should be interpreted with caution for these three GRBs. 

(iii) GRBs 1025 and 8086, both long duration bursts, have been found to have very low 
variabilities and, therefore, low luminosities. This results in very high Lorentz factors 
in our calculation (20,000 and 5,000, respectively). Given the uncertainties involved 
in the variability calculation and the variability-luminosity relationship, we feel that 
these luminosities and Lorentz factors are unlikely to accurately represent these bursts. 
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7. Results and Discussion 

As already mentioned, our results are summarized in a series of tables. They summarize 
our fits to the BATSE data that provide the input spectrum for calculating the neutrino 
emission. The neutrino event rates for the 4 classes and the 4 models are also tabulated (see 
tables 1 through 8). In tables 9-12, we summarize the event rates for a generic kilometer- 
scale telescope, such as IceCube, for the 4 classes of GRBs. For the duration and in the 
direction of a GRB, the background in a neutrino telescope should be negligible. Therefore 
individual events represent a meaningful observation when coincident with a GRB. We next 
evaluate the prospects for such observations. 

The AMANDA experiment has operated for approximately four years (1997-2000) that 
overlap with the BATSE mission (Barouch & Hardtke 2001). Data has been collected for 
several hundred GRBs with an effective telescope area of order 5,000 m 2 . AMANDA-II, the 
completed version of the experiment, with approximately 50,000 m 2 effective area for the 
high energy neutrinos emitted by GRBs, was commissioned less than half a year before the 
end of the BATSE mission. Nevertheless, of order 100 bursts occurred during that period 
for which coincident observations were made. With effective areas significantly below the 
reference square kilometer of future neutrino telescopes, AMANDA is not expected to test 
the GRB models considered here. For long duration GRBs, the most common classification, 
we anticipate on the order of 0.01 events (muons+showers) per square kilometer per GRB 
for models 1 and 2. For AMANDA-II, with one twentieth of this area, and only capable 
of observing northern hemisphere GRBs (~ 500/yr), we predict on the order of .3 events 
(muons+showers) per year of observation. AMANDA-B10, with smaller area, should observe 
one tenth of this rate. It is interesting, however, to note that such experiments are on the 
threshold of observation at this time. IceCube with a square kilometer of effective area, now 
under construction, will likely cross this threshold. 

For the high energies considered here, IceCube should be able to make observations of 
1000 bursts over 4n steradians during one year. Where models 1 and 2 are concerned, we 
expect that an event (muon or shower) will be observed from roughly 10 GRBs. For models 
3 and 4 we predict ~ 7 events per year for models with t s d ~ 0.07 yr (model 4) , but only 
around one event per ten years if t S( j is somewhat larger, such as 0.4 yr (model 3). The 
prospects for observation depends strongly on t sd , as expected. 

We expect that classes of GRBs with lower average fluence, such as short duration 
GRBs and X-ray flash candidates, will be more difficult to observe. X-ray flash candidates, 
although unlikely to be observable in models 1 and 2, could possibly be observed in models 
3 and 4. If t sd ~ 0.07 yr (model 4), we predict order one event (muon or shower) for every 
1000 bursts. With only ~ 100 X-ray flashes thought to occur per year, such an observation 
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will be unlikely. If such events are more common, however, they may contribute significantly 
to the diffuse high energy neutrino flux. Observations of neutrinos from low peak energy 
GRBs would be strong evidence for a supranova progenitor model. 

It is interesting to note that the majority of the neutrino events from GRBs come from 
a relatively small fraction of the GRB population (Alvarez, Halzen & Hooper 2000; Halzen 
& Hooper 1999). In Figure 1 we show some of the factors which go into this conclusion. 
The occasional nearby (low redshift), large /„., high fluence and/or near horizontal burst can 
dominate the event rate calculation. The distribution of the number of events per burst (or 
X-ray flash candidate) is shown in Figure 2. 

In summary, taking advantage of the large body of GRB statistics available in the 
BATSE catalog, we have attempted to estimate the neutrino fluxes and event rates in neu- 
trino telescopes associated with GRBs for a variety of theoretical models. Our analysis has 
yielded several conclusions: 

(i) Gamma-ray bursts with high fluence, most often long duration bursts, provide the best 
opportunity for neutrino observations. 

(ii) For typical gamma-ray bursts, proton interactions with fireball photons provides the 
largest neutrino signal. We have also shown that the rates are relatively model inde- 
pendent. 

(iii) For gamma-ray bursts with very low peak energies, possibly associated with X-ray 
flashes, very little energy is transfered into pions (and, therefore, neutrinos) in inter- 
actions with fireball photons. Interactions with a surrounding pulsar wind bubble, 
however, can yield interesting neutrino fluxes. This is an illustration that observation 
of neutrinos is likely to help decipher the progenitor mechanism. 

(iv) While our calculations indicate that existing neutrino telescopes, such as AMANDA, 
are not likely to have the sensitivity to observe gamma-ray burst neutrinos, next gen- 
eration, kilometer-scale observatories, such as IceCube, will be capable of observing on 
the order of ten bursts each year. 
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A. Appendix: Photomeson Interactions of Protons With GRB Photons 

In this appendix, we describe the production of neutrinos in interactions of protons 
and photons in the GRB fireball. Protons predominantly produce the parent pions via the 

processes 

— > A -> mr + (Al) 

and 

jrf -> A -> pvr (A2) 

which have very large cross sections of a a ~ 5 x 10~ 28 cm 2 . The charged 7r's subsequently 
decay producing charged leptons and neutrinos, while neutral 7r's decay into high-energy pho- 
tons. For the center-of-mass energy of a proton-photon interaction to exceed the threshold 
energy for producing the A-resonance, the comoving proton energy must meet the condition: 

2 2 

mi — mi 

4 > ^r^- (A3) 



7 



Throughout this paper, primed quantities are measured in the comoving frame and unprimed 
quantities in the observer frame. In the observer's frame, 



r 2 

e p > 1.4 x 10 16 ^^eV, (A4) 

^7,MeV 



resulting in a neutrino energy 



e u = -Ax^ )e p >7x 10 14 -^eV, (A5) 

4 £ 7 ,MeV 

where I\5 = T/10 2 ' 5 is the plasma expansion (bulk) Lorentz factor and £ 7i mcV — £ 7 /lMeV is 
the photon energy. (x p ^ n ) ~ 0.2 is the average fraction of energy transferred from the initial 
proton to the produced pion. The factor of 1/4 is based on the estimate that the 4 final state 
leptons in the decay chain n + — > — > z/^e+z/gZ/^ equally share the pion energy. These 
approximations are adequate given the uncertainties in the astrophysics of the problem. 

For each proton energy, the resulting neutrino spectrum traces the broken power law 
spectrum of photons which we fit to the BATSE data using the broken power law parame- 
terization 



(A6) 
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Summing over proton energies results in a neutrino spectrum with the same spectra slopes, 
a and f3, as for the gamma-ray spectra in the BATSE data, but with a break energy of order 
1 PeV in the observer frame: 

4 = 7x10" 1 J™ eV . (A7) 

We here explicitly introduce the dependence on source redshift, z. The highest energy 
pions may lose some energy via synchrotron emission before decaying, thus reducing the 
energy of the decay neutrinos. The effect becomes important when the pion lifetime ~ 
2.6 x 10~ 8 e' n I \m n c 2 ) s becomes comparable to the synchrotron loss time 

^ = Aarmle^ (A8) 

where U' B = B' 2 /8tt is the energy density of the magnetic field in the shocked fluid. e#, 
the fraction of the internal energy carried by the magnetic field, is defined by the relation 
AnR 2 cY 2 B' 2 /87r = e^-Lint, where R ~ 2T 2 ct v , is the collision radius. The collision radius R 
is obtained from the consideration that different shells in the shocked fireball have velocities 
differing by Av ~ c/2r 2 , where T is an average value representative of the entire fireball. 
Different shells emitted at times differing by t v therefore collide with each other after a time 
t c ~ ct v / Av , i.e. at a radius R = ct c ~ 2T 2 ct v . A detailed account of the kinematics can be 
found in Halzen & Hooper 2002. We can now compare the synchrotron loss time with the 
time over which the pions decay: 

t' sy J4 = 0.21 e e e^L^ri 5 tl_ 2 e- 2 8 . (A9) 

Here e e is the fraction of internal energy converted to electrons, t v -2 = t v /10~ 2 s is the time 
scale of fluctuations in the GRB lightcurve, L l 52 — L^/10 52 erg/s is the 7-ray luminosity of 
the GRB and e n ±s = 5,r/10 18 eV, is the pion energy. In deriving Eq.(A9) we have assumed 
that the wind luminosity carried by internal plasma energy, L int , is related to the observed 
7-ray luminosity through L int = L 7 /e e . This assumption is justified because the electron 
synchrotron cooling time is short compared to the wind expansion time and hence electrons 
lose all their energy radiatively. 

The radiative losses become important for t' sync < r^, which corresponds to e n > e ws ~ 
4e us , where 

^ M = ^ ^B /2 L-% 2 Tt,t^ 2 eV. (A10) 

Neutrinos from muon decay have a lifetime 100 times longer than pions, the energy cutoff 
will therefore be 10 times smaller: 

1 D 16 

*U = TT~z e e /2e B /2L ^i. 5 t v ,-2 eV. (All) 
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Above this energy, the slope of the neutrino spectrum steepens by two to (f3 + 2). 

To normalize the neutrino spectrum to the observed GRB luminosity, we must calculate 
the fraction, of fireball proton energy lost to pion production. The fraction of energy 
converted to pions is estimated from the ratio of the size of the shock, AR', and the mean 
free path of a proton for photomeson interactions: 

AR' 

f„~—(x p ^). (A12) 

Here, the proton mean free path is given by A P7 = l/n 7 <7A, where n 1 is the number density 
of photons. The photon number density is given by the ratio of the photon energy density 
and the photon energy in the co moving frame: 

< A13 > 

Using these equations, and recalling that R ~ 2T 2 ct v , we obtain that 

^7 \ /(<h\ = ^0 

t v T 5 AR' ) I \ V ) 16ncH v r ^ , , . 

and the fraction of proton energy converted to 7r's is 



Ul ~ ( 167rc 2 t„r 5 AR' ) / ( T ) 16ncH„T 4 AR'e^ ' (M4) 



h- p 4 , lfi 2 ~U.2X . (Albj 

fc 7 1 t„ 1U/IC 1 2.5t t)) _2fc 7 ]yj e v 



This derivation was performed for protons at the break energy. In general, 

Ue P ) - 0.2 L ^ x { *>J , (A16) 

1 2.5 t f,-2 fc 7 ,McV L \ t P/ t p) t P ^ t p 

where e b p is given by Eq.(A4). 

As we can see from Eq.(Al5), f n strongly depends on the bulk Lorentz factor T. It has 
been pointed out by Halzen & Hooper (1999) and Alvarez, Halzen & Hooper (2000) that, 
if the Lorentz factor T varies significantly between bursts, then the resulting neutrino flux 
will be dominated by a few bright bursts with f n close to unity. However, Guetta, Spada & 
Waxman (2001a) have shown that burst-to-burst variations in the fraction of fireball energy 
converted to neutrinos are constrained. First, the observational constraints imposed by 7-ray 
observations, in particular the requirement e b > 0.1 MeV, imply that wind model parameters 
(r, L, t v ) are correlated (Guetta, Spada & Waxman 2001b) and that T is restricted to values 
in a range much narrower than AT/T ~ 1. For instance, for values of 7 much smaller than 
average the fireball becomes very dense with abundant neutrino production. Such fireballs 
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will also produce a thermal photon spectrum which is not the case for the events considered 
here. Second, for wind parameters that yield f n values significantly exceeding 20%, only a 
small fraction of pion energy is converted to neutrinos because of pion and muon synchrotron 
losses as can be seen from E q.(AlO). 

We will use two methods to determine the value of the bulk Lorentz factor, T. For bursts 
with high break energies, > 500 keV, T cannot differ significantly from the minimum value 
for which the fireball pair production optical depth is ~ 1 near the maximum energy of 7- 
rays produced, e 1 ^ max . EGRET has observed 7-rays with energies in excess of 1 GeV for six 
bursts, although the maximum 7-ray energy should be lower for the majority of GRBs. We 
choose 100 MeV as the default value, therefore, 



250 



Note that in the end, the value of the Lorentz factor depends weakly on luminosity, time 
structure and maximum 7-ray energy. 

For bursts with lower break energies Eq.(A17) may not be reliable because the Lorentz 
factor of these GRBs may be larger than estimated. Guetta Spada & Waxman (2001b) have 
argued that the X-ray flashes identified by BeppoSAX could be produced by relativistic 
winds where the Lorentz factor is larger than the minimum value given in Eq.(Al7) required 
to produce a GRB with the characteristic photon spectrum. For GRBs with low break 
energy we, instead, relate the Lorentz factor to the peak energy of the 7-ray spectrum. The 
characteristic frequency of synchrotron emission is determined by the minimum electron 
Lorentz factor 7 m pa e e {m p /m e ) and by the strength of the magnetic field given above, 
before Eq.(A9). The characteristic energy of synchrotron photons, = Yhj^eB' /2nm e c, at 
the source redshift is 

T -1/2 

4 - i^i^f-MeV. (A18) 

1 2.5 l v -2 

For bursts with e b ^ < 500 keV we will evaluate T from the break photon energy given above. 
At present no theory allows the determination of the values of the equipartition fractions e e 
and 6b- Eq.(Al8) implies that fractions not far below unity are required to account for the 
observed 7 ray emission and this is confirmed also by simulations (Guetta Spada & Waxman 
2001b). For bursts with very large values of T, the peak energy is shifted to values lower than 
40 keV. This could be the X-ray bursts detected by BeppoSAX (Guetta Spada & Waxman 
2001b). From Eq.(Al5), we estimate that the neutrino flux from such events is expected to 
be small. 



We have now collected all the information to derive the neutrino from the observed 
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7-ray fluency F 7 : 

if/" ~ 8^Rio) A - (A19) 

BATSE detectors measure the GRB fluence F 7 over two decades of photon energies, ~ 0.02 
MeV to ~ 2 MeV, corresponding to a decade of energy of the radiating electrons. The factor 
1/8 takes into account that charged and neutral pions are produced with roughly equal 
probabilities, and each neutrino carries ~ 1/4 of the pion energy. Using Eq.(Al6), 

i I p r f {e v /e*)f> e v < e\ 

eldN v /de v « 0.2 ^ x ( e „/4)« 4 <e v <el , (A20) 

8 e e ln(10) r 16 ^ v { ^^^-i ^ > £ , 

where £^ and are given by Eq.(A7), Eq.(AlO) and Eq.(All). This spectrum is shown in 
Fig.(3). 

This result depends on a number of somewhat tenuous assumptions. Simulations (GSW) 
actually suggest that one can simply fix f w ~ 0.2 at the break energy and derive the v flux 
directly from the 7-flux. Doing better may require a better understanding of the fireball 
phenomenology than we have now. First, the variability time may be shorter than what is 
observed; in most cases variability is only measured to the smallest time scale that can be 
detected with adequate statistics. Second, the parameters e e and e# are uncertain. Third, 
the luminosity-variability relation used to derive the luminosity for bursts with no measured 
redshift (see §4.1) is uncertain, and in addition may have large fluctuations around the pre- 
diction. We have therefore decided to do the detailed analysis described above as well as 
an alternative analysis that assumes f n = 0.2, at the break energy, for all bursts and deter- 
mines the neutrino flux directly from the observed gamma-ray fluence. For this alternative 
approach, 

rlN fl 9 TP { ( £l/ / £ k) E v < e u 

^Z, ,,^ 4<e : <K ' (A21) 

K J { (e u /e b u ) a (e u /eir 2 e u > e s u 

We refer to the models based on Eq.(A21) and Eq.(A20) as models 1 and 2, respectively. In 
Fig. 3 we show the muon neutrino spectrum for our fiducial parameters in models 1 and 2. 



B. Appendix: Photomeson Interactions of Protons with External Photons in 

the Supranova Model 

In this section, we consider the neutrino photoproduction on external photons in supra- 
nova GRBs (Guetta & Granot 2002a). The external radiation field surrounding the final 
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black hole is referred to as the pulsar wind bubble (PWB). A pulsar wind bubble is formed 
when the relativistic wind (consisting of relativistic particles and magnetic fields) that em- 
anates from a pulsar is abruptly decelerated (typically, to a Newtonian velocity) in a strong 
relativistic shock, due to interaction with the ambient medium. 

Fractions of the post-shock energy density go to the magnetic field, the electrons and 
the protons, respectively. The electrons will lose energy through synchrotron emission and 
inverse- Compton (IC) scattering. In (Guetta & Granot 2003) a deep study of the character- 
istic features of the plerion emission has been carried out. As shown in (Guetta & Granot 
2003), the electrons are in the fast cooling regime for relevant values of t s &, and therefore most 
of the emission takes place within a small radial interval just behind the wind termination 
shock. 

The mechanism for neutrino production through photomeson interaction with the ex- 
ternal photons dominates when the lifetime of the supramassive neutron star, t s d < 0.2 yr 
(or t s d < 2 yr for X-ray flashes). Neutrinos generated in this way have typical energies 
e u ~ 10 15 - 10 17 (10 19 ) eV (10 19 eV for X-ray flashes). As in models 1 and 2, these neutri- 
nos are emitted simultaneously with the prompt 7-ray (X-ray) emission. For even shorter 
lifetimes t sd < 0.1 yr, the z/'s would not be accompanied by a detectable GRB because the 
Thomson optical depth on the PWB is larger than unity. 

As before, protons of energy e p interact mostly with photons that satisfy the A-resonance 
condition, e p ^ = 0.3GeV 2 /e 7 , where, in this the PWB photon energy. The 

minimum photon energy for photomeson interactions corresponds to the maximum proton 
energy £ 7 , m i n = 0.3 GeV 2 /e Pimax ~ 3 x 10~ 3 (£ Pima x/T0 20 eV) _1 eV. For reasonable model 
parameters (t S( j > 0.026 yr) this energy exceeds self absorption frequency of the PWB 
spectrum (Guetta & Granot 2003). Moreover, for t S( j < 12 yr the electrons are in fast cooling 
and emit synchrotron radiation. Therefore the relevant part of the spectrum consists of two 
power laws, <in 7 /cfe 7 oc e 7 3//2 for e 1 < = hvbm and <in 7 /cfe 7 oc e 7 S//2 1 for e 1 > £ 7 &, where 
Vbm ~ 1-6 x 10 15 t^ d 3 ( 2 1 Hz is the peak frequency of the PWB spectral energy distribution 
(z/F^), n 7 is the number density of photons and s m 2.2 is the power law index of the PWB 
electrons. 

The normalization factor of the target photon number density is determined by equating 
the pulsar wind luminosity in pairs, £ e -Erot/4d (£e is the fraction of the pulsar wind energy 
in the e ± component and E TOt is the total energy of the pulsar wind) to the total energy 
output in photons, which is ~ U ph A^R 2 at R 3> R s = fRb, where R s is the radius of 
the pulsar wind termination shock, which is a factor, /, smaller than the PWB radius, 
Rb- At R < R s , which is relevant for our case, the photons are roughly isotropic and 
U p h becomes roughly constant, and assumes the value U p h ~ (, e E ro t /4d2vr(/ 'Rb) 2 c. As we 
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mentioned above, the relevant target photons for photomeson interaction with high energy 
protons are the synchrotron photons and the fraction of the total photon energy that goes 
into the synchrotron component is U syn /U p h = (1 + Y b ) y 'e\>B / 'e&e, where e?, e and e&£ are 
the fractions of the PWB energy in the electrons and in magnetic field, respectively. 

This spectrum is consistent with the spectrum of known plerions like the Crab, that can 
be well fit by emission from a power law distribution with a power law index value ~ 2.2. 
In our case, the only difference is that there is a fast cooling synchrotron spectrum, rather 
than a slow cooling one. There are no observations of the spectrum from very young plerions 
where there is a fast cooling spectrum (as they are more rare). 

This radiation will be typically hard to detect (for t s d < 1 yr and z ~ 1), but might be 
detected for closer (though rarer) PWBs. 

The proton energy satisfying the A-resonance condition with PWB photons of energy 

p 2 n3/2 .3/2 

^ = 4-4x10 eV, (Bl) 

'/2/3 e 6e/3 e feB, -3-^53 7to,4.5 

where t s d,-i = ^sd/0.1 yr, = /3 b /0.1 is the velocity of the SNR shell (in units of c), 
E 53 = E rot /10 53 erg, 7^4.5 = 7 W /10 4 - 5 is the Lorentz factor of the pulsar wind, £ e = £ e /3/3, 
e^e = ebe/3/3, €bB = 10~ 3 e5£,-3) and 77 = (2/3)772/3 is the fraction of the wind energy that 
remains in the PWB. The latter is the fraction that goes into the proton component and 
is, unlike the electron component, not radiated away. The corresponding neutrino energy is 
£1/6 ~ £pf>/20 ~ 2 x 10 15 t^ 2 _ 1 eV. For i S( j ~ 0.1 yr, this energy is similar to those obtained in 
interactions with GRB photons in the previous section. Photon emission is only detected in 
coincidence with these neutrinos if the Thompson optical depth is < 1, which is the case for 
t sd ~ 0.1 yr and a clumpy SNR (Guetta & Granot 2003; Inoue, Guetta & Pacini 2003). In 
the case of a uniform shell the condition is t sd > 0.4 yr, corresponding to e vh > 2 x 10 16 eV. 

As before, the internal shocks occur over a distance R = 2T 2 ct v . Thus, the optical depth 
for photo-pion production by protons of energy e p , is 

drij _ l-0£f/3#53 2 r|.5^ -liSp/epb) 13 

Tjn - < Vy £ 7^~- K - f2 ^5/2^5/2 ^ 2 ^1/2 ^3/2 ^ B2 J 



^l/3'/2/3 e fee/3 IpA-^Pb -l'-sd-l 



where r 2 . 5 = T/10 2 - 5 , t v - 2 = K/10' 2 s, a P1 « 5 x 10~ 28 cm 2 , e 7 = 0.3 GeV 2 /5 p , / = / V3 /3 
and P is the spectral slope of the seed PWB synchrotron photons with (3 = s/2 (1/2) for 
£ P < £ P b (z P > £pb)- The fraction of the proton energy that is lost to pion production is given 
by 

fpn(£ P ) ~ 1 - exp [-T P7 (£ p )/5] « min [1, r P7 (£ p )/5] . (B3) 
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The factor of 5, as before, takes into account that the proton loses ~ 0.2 of its energy 
in a single interaction. We denote e p for which / p7r (e p ) ~ 1 by £ pr i8 = £ pr /10 18 eV [i.e. 
r P7 (e pr ) = 5], and obtain 



n onf 20 / 11 25/11 3/11 „43/22 63/22 
U - ZU Jl/3 ^2/3 £ i>e/3 P h,-l sd,-l £■ 



,8/11 1/2 „21/22 2/11 r 20/ll 10/11 £ 
?e/3 e bB, -3^53 'p.4.5 1 2.5 r v,-2 p0 

1 2f 4 n 5/ V -v 2 fl 5/2 t 9/2 

1 ' Z Jl/3 f '2/3 i; i)e/3 'p,4.5Pb,-l r sd,-l £ pT ^ 1 
^e/3 e f)B,-3-' 1 '53 1 2.5 t i>,-2 p0 
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As before, the decay of charged pions created in interactions between PWB photons and 
GRB protons, produces high energy neutrinos, 7r + — > ji + + z/^ — > e + + z/ e + z/^ + z/ M , where 
each neutrino receives ~ 5% of the proton energy. 

The total energy of the protons accelerated in the internal shocks is expected to be 
similar to the 7-ray energy produced in the GRB (Waxman 1995). This implies a fluence, 

U = fof P u , /o = |^ = lxlO- 5 §^^|, (B5) 
M 32nd z L d^ L28 cm 2 

_ / de p (dN p /de p )e p f pu (e p ) 
pv J de p {dN p /de p )e p 

where £ 7i53 = £ 7jiso /10 53 erg is the isotropic equivalent energy in 7-rays, f pv {e p ) = f p7 ,{e p )f^ u (e p ), 
while fp-n^Sp) is given in Eq.(B3) and f^Sp) is the fraction of the original pion energy, 
En ~ 0.2e p , that remains after decay. 

The pions may lose energy via synchrotron or inverse Compton (IC) emission. If these 
energy losses are significant, then the energy of the neutrinos will be reduced as well. Follow- 
ing arguments already presented in appendix A, we find that /^(^p) ~ 1 — exp(— t' rad /T^) « 
min(l, t f rSud /rl), where t' % ~ 2.6 x 10~ 8 e' n / (m n c 2 ) s is the lifetime of the pion, and t' rSud = 
(t'~n + ^c 1 ) -1 ~ mm (^syn^ic) * s the time for radiative losses due to both synchrotron and 
IC losses. The time, t' syD , is given by Eq.(A8) and 

t{c = tojmlelufa) ' (B6) 

where U^e^) is the energy density of photons below the Klein-Nishina limit, e ly KN = 
(m^c 2 ) 2 1 'e w . IC losses due to scattering of the GRB photons were shown to be unimportant 
(Waxman & Bahcall 1997). We therefore only consider the IC losses from the upscattering 
of the external PWB photons, and find 
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where e e and e# are the equipartition parameters of the GRB, and e^is = £tt/10 18 eV. 
The radiative losses become important for t' Tad < t£, which corresponds to e p > e ps = 
mm(e p y s n ,s p c s ) w 5e ns « 20e„ s , where 



The protons may also lose energy via p — 7 interactions with the GRB photons (Waxman 
& Bahcall 1997). However r P7 for this process is typically < 1, so that it does not have a 
large effect on yry interactions with the PWB photons, on which we focus. 

Because the lifetime of the muons is ~ 100 times longer than that of the pions, they 
experience significant radiative losses at an energy of e^ s ~ e ns /10 ~ e ps /50. This causes a 
reduction of up to a factor of 3 in the total neutrino flux in the range ~ (0.1 — l)e us , since 
only that are produced directly in n + decay contribute significantly to the neutrino flux. 
Note that since both ratios in Eqs. (A9) and (B7) scale as e~ 2 , we always have 4ad/ r ^ ^ 6 7^ 
and therefore, the spectrum steepens by a factor of (e v /s vs )~ 2 for e v > e vs . This is evident 
because f wu (e p « 5e w ) ss min(l, t' sy Jr^, t' ic /r^). 

In Figure 4 we show the proton energies that correspond to the neutrino break energies 
£i/sj e ub and w £ pr /20, as a function of t sd . From Eq. (B2), (B8) and (B9), we conclude 
that t p1 oc T 2 ^ 3 ^ 2 , while e s p y s D oc T 4 and e 1 ^ oc t% 2 . For a fixed value of t s d -1 = 1, ^.5 > 1 
implies r P7 > 1 and e ps = e l p c s = constant ~ 10 18 eV, while r 2 .5 < 1 implies e ps = e s p y s n oc T 4 . 
This implies increased neutrino emission reaching higher energies for larger values of T. Since 
large T implies lower synchrotron frequency for the prompt GRB, this may be relevant to 
X-ray flashes, assuming that they are indeed GRBs with relatively large Lorentz factors 
and/or a large variability time, t v (Guetta, Spada & Waxman 2001b). 

As can be seen from Figure 4, depending on the relevant parameters, there are four 
different orderings of these break energies: (i) e VT < e u b < £us, (ii) &vb < £vt < £us, (hi) 
£1/6 < £vs < £vt and (iv) e vs < e v b < e UT . Each of these results in a different spectrum 
consisting of 3 or 4 power laws. Analytically these spectra are: 




(B8) 
(B9) 




(B10) 
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.2 dNy 
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ub \ 2 ( Svt\ ' 

UT ) \ £v ) 
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For our analysis we consider two characteristic values of t s a, 0.07 yr and 0.4 yr, and refer to 
the models as 3 and 4, respectively. In the case of Model 3, the GRB is seen only if the shell 
is sufficiently clumpy while in model 4 the GRB should always be detectable. 

The muon neutrino spectrum is shown in Figure 5 for our fiducial parameters and 
tsd = 0.01, 0.07, 0.4, 30 yr. The spectrum of the other neutrino flavors is the same. 



C. Appendix: Calculations of Event Rates 

A compilation of the probability that a GRB neutrino is actually detected as a muon, a 
tau or a shower by an underground detector is shown in Fig. 4 as a function of the neutrino 
energy. These are required to convert neutrino spectra from GRBs to event rates. We present 
in this appendix the formalism for doing this conversion. 
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C.l. Showers 

The number of shower events in an underground detector from a neutrino flux 
produced by a single GRB with duration T is given by 

. r w<t> 

N sh = Y. 2nAT J dE Uil ^(E Ui )P smv (E Vi ,e z ) 



1,3 

i.i 

1 da J 
'aJ(E~j^y~ 



x/ dy—ypr^ ~ (E Vi , y) P mt (E Ui , y,0 z ) , (CI) 



*min 



where 9 Z is the zenith angle (9 Z = is vertically downward). The sum is over neutrino (and 
anti-neutrino) flavors i — e, p,T and interactions j = CC (charged current) and NC (neutral 
current). A is the detector's cross sectional area with respect to the v flux, and d$> u JdE v . 
is the differential neutrino flux that reaches the Earth. For % — r, CI is modified to include 
the effects of regeneration of neutrinos propagating through the Earth, as will be discussed 
further on. 

P surv is the probability that a neutrino reaches the detector, i.e. is not absorbed by the 
Earth. It is given by 

P surv = exp[-X(^) ( x tot (^)iV A ] , (C2) 
where Na — 6.022 x I0 23 g _1 , and the total neutrino interaction cross section is 

(T tot = (T CC + C7 NC (C3) 

This is somewhat conservative because it neglects the possibility of a neutrino interacting via 
a NC interaction and subsequently creating a shower in the detector. X(6 Z ) is the column 
density of material a neutrino with zenith angle 6 Z must traverse to reach the detector. It 
depends on the depth of the detector and is given by 

X(0z)= [ P (r(9 z ,l))dl , (C4) 

the path length along direction 9 Z weighted by the Earth's density p at distance r from the 
Earth's center. For the Earth's density profile we adopt the piecewise continuous density 
function p(r) of the Preliminary Earth Model (Dziewonski 1989). 



Pint is the probability that the neutrino interacts in the detector. It is given by 

L 



P int = 1 - exp 



(C5) 



where, for showers, L is the linear dimension of the detector, and D(E Ui ) is the mean 
free path for neutrino interaction of type j. For realistic detectors, L <^ L?{E U% ), and so 
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P mt L/LP{E Vi ). To an excellent approximation the event rate scales linearly with detector 
volume V = AL. 

The inelasticity parameter y is the fraction of the initial neutrino energy carried by the 
hadronic shower (rather than the primary lepton). The limits of integration depend on the 
type of interaction and on the neutrino flavor. For NC u e interactions and all and u T 
interactions, y max = 1 and y min = E^ r /E u , where E^ is the threshold energy for shower 
detection. For CC u e interactions, the outgoing electron also showers, therefore y max = 1 
and y min = 0. 



C.2. Muons 



Energetic muons are produced in CC interactions. For a muon to be detected, it 
must reach the detector with an energy above its threshold E^ T . The expression of Eq.(Cl) 
then also describes the number of muon events after the replacement 



int 



cxp 



RfijE^, 9 Z ) 
L CC (E^) 



(C6) 



where is the range of a muon with initial energy 



y)E Vli and final energy E 1 ^ 1 . 



We will assume that muons lose energy continuously according to 

dE 

7x = - a -' 3E - 



(C7) 



where a = 2.0 MeVcm 2 /g and j3 = 4.2 x cm 2 /g (Dutta et al. 2001). The muon range 
is then 

a + PE^ 



a + (3E^ 



(C8) 



In this case, y max = 1 - E^ T /E„ and y min = 0. 



The event rate for muons is enhanced by the possibility that muons reach the detector, 
even if produced in neutrino interactions kilometers from its location. Note, however, that 
this enhancement (i.e. is # 2 -dependent: for nearly vertical down-going paths, the path 
length of the muon is limited by the amount of matter above the detector, not by the muon's 
range. This is taken into account in the simulations. Fig. 7 shows the probability of detecting 
a muon generated by a muon neutrino as a function of the incidence zenith angle. The figure 
illustrates the effect of the limited amount of matter above the detector, as well as the 
neutrino absorption in the Earth. Absorption is not important at E v — 1 TeV, and the 
muon range at this energy (for a muon energy threshold of 500 GeV as we adopted in the 
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figure) is not limited by the amount of matter above the detector. As a consequence the 
probability at 1 TeV is weakly dependent on zenith angle. Absorption in the Earth starts 
to be important for upgoing neutrinos of energy above ~ 100 TeV - 1 PeV, and restricts the 
neutrino observation to the horizontal and downward directions at extremely high energy 
(EeV range) as can be seen in the figure. The large muon range at PeV and EeV energies, 
much larger than the depth at which the detector is located (1.8 km vertical depth), reduces 
considerably the detection probability of downgoing neutrinos. 

The largest background to a GRB signal consists of muons from atmospheric neutrinos. 
However, for GRB observations, the time and angular windows are very small and this 
background can be easily controlled, as we will illustrate. Following (Dermer & Atoyan 
2003), the number of background events is approximately given by 



iV, 



b g 



J cM J dtj™ dE v ^-P v ^{E v ) (C9) 



where Q and t are the solid angle and time considered, respectively, and P U ^^(E U ) is the 
probability of a muon neutrino generating a muon in the detector volume. IceCube is 
designed to have angular resolution smaller than 1 degree at the relevant energies. We 
consider a 1 degree cone for the solid angle calculation. Considering a long burst of duration 
~ 100 seconds, and using the following approximate atmospheric neutrino spectrum 

dN u f 5 x lO-^^GeVcm-^sr- 1 E v > 10 5 GeV 

(C10) 



dE v \ 1.6 x lO-^-^GeVcm-^sr- 1 E v < 10 5 GeV ' 

we find that for a natural threshold (minimum energy) of ~ 100 GeV, we expect ~ .003 
background events per burst. More practically, an energy threshold in the range of 1-10 TeV 
could be imposed which would reduce this background by an addition factor of 50 to 2500, 
respectively. For a naive illustration, consider one years of observation, with 1000 bursts, 
each of duration of 10 seconds and a 1 TeV energy threshold imposed. For such a example, 
less than 0.01 total background events are predicted. 



C.3. Taus 

Taus are produced only by CC v T interactions. This process differs significantly from 
the muon case because tau neutrinos are regenerated by the production and subsequent tau 
decay through v T — > r — > v T (Halzen & Saltzberg 1998). As a result, for tau neutrinos, CC 
and NC interactions in the Earth do not deplete the v T flux, they only reduce the neutrino 
energy down to a value where, eventually, the Earth becomes transparent. We implement 
this important effect using a dedicated simulation that determines E Vt (E Vt ,9 z ), the average 



energy of the v T when reaching the detector. It depends on the initial energy E Vt and zenith 
angle 9 Z . The tau event rate is then given by 



N T = 2.AT J d E„ T § 
xP int e((l-y)E 



{E VT ) / dy 



/*ymax 



1 



da cc _ 



(Eu T ,y) 




° cc {E Vt ) dy 



(Cll) 



P int depends on the geometry of the neutrino tau induced event. For events consisting on a 
minimum-ionizing track going through the detector it is given by: 



where R T ((1 — y)E Ur ) is the range of the produced tau evaluated at the energy of the tau 
after regeneration. R T is given by Eq.(C8) with (5 = 3.6 x 1CT 7 cm~ 2 /g (Dutta et al. 2001). 
The last factor takes into account the requirement that the tau track be long enough to be 
identified in the detector. We require P* hr ~ 2.5 x 10 6 GeV so that the tau decay length 
is larger than the 125 m string spacing in IceCube. It is not clear with what efficiency 
through-going tau events can be separated from low energy muons. 

Those tau events that include one shower (lollipop events) or two showers (double bang 
events) inside the detector volume will be identifiable. For these cases we use P int obtained 
by a dedicated simulation that determines the probabilities for double bang and lollipop 
geometries shown in Fig. (6). The rate of downgoing lollipop events in a km 3 neutrino tele- 
scope is expected to be of the order of the rate of down-going shower events, probably slightly 
smaller. Double bang events will be mostly observed for neutrino energies in a limited range 
between roughly 10 and 100 PeV. 

Finally, those events in which the tau decays into muons that reach the detector are 
counted as muon events. Pi nt is given by Eq.(C12) where we use as range the sum of the 
range of the tau evaluated at the regeneration energy and the range of the muon at the 
energy that carries in the decay of the tau. 

As with muons, at very high energies taus can travel several kilometers before decaying 
or suffering significant energy loss. The enhancement to tau event rates from this effect is 
^-dependent as discussed above for muons. 




(C12) 
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Table 1: Characteristics for GRBs with redshifts measured from optical afterglow observa- 
tions. The fluence and zenith angle of each burst are taken from the BATSE catalog, whereas 
the luminosity is derived from the redshift and fluence. The break energy is obtained by a fit 
to the BATSE data and the Lorentz factor is calculated as described in Appendix A. Note 
that GRB 6707 is an anomalous burst (see §6). 



BATSE # 


F 7 (erg/cm 2 ) 


L 7 (erg/s) 


z 


£$(MeV) 


r 


9 Z (degrees) 


8079 


1.61 10" 


-6 


4.9 10 51 


1.118 


0.460 


222. 


166.2 


7906 


2.51 10" 


-4 


2.0 10 52 


1.020 


1.460 


281. 


100.9 


7648 


5.83 10- 


-6 


2.0 10 49 


0.434 


0.710 


89. 


101.8 


7560 


2.06 10- 


-5 


2.0 10 51 


1.619 


0.100 


925. 


9.8 


7549 


2.11 10" 


-4 


1.2 10 51 


1.300 


1 


176. 


63.8 


7343 


4.87 10" 


-4 


4.0 10 51 


1.600 


0.280 


214. 


132.3 


6891 


6.22 10- 


-5 


1.5 10 49 


0.966 


0.280 


85. 


102.0 


6707 


4.01 10" 


-6 


1.3 10 45 


0.009 


0.430 


25299. 


36.9 


6533 


1.25 10" 


-5 


2.0 10 52 


3.420 


0.460 


281. 


156.2 


6225 


3.96 10- 


-6 


1.2 10 50 


0.835 


1.400 


119. 


170.6 
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Table 2: Estimated neutrino event rates for GRBs with redshifts measured from optical 
afterglow observations. Rates for muon and shower events are shown for the IceCube neutrino 
telescope. Tau events are not shown. Note that GRB 6707 is an anomalous burst. Model 

1 and model 2 correspond to the internal shock scenario. In model 1 we fix the fraction 
of energy transferred from protons to pions to a constant value f w = 20%, and in model 

2 we calculate f w from the observed GRB characteristics using Eq.(Al6). Models 3 and 4 

correspond to the supranova scenario, with choices of the relevant parameter t s< j = 0.4 and 

0.07 years respectively. Note that in the supranova scenario (models 3 or 4), internal shock 

processes (models 1 or 2) also contribute to the event rate. The total predicted rate for 

models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the rates 

shown for the internal shock model. 

Model 1 Model 2 Model 3 Model 4 



BATSE # 




n 




Shower 


n 




Shower 


n 




Shower 




n 




Shower 


8079 


3.1 


10" 


-5 


6.8 10~ 6 


5.2 10" 


-5 


1.2 10~ 5 


1.7 10- 


-8 


2.9 10~ 9 


3.1 


10" 


-6 


5.8 10" 7 


7906 


9.9 


10" 


-2 


2.1 10~ 2 


9.7 10" 


-2 


2.1 10- 2 


9.1 10- 


-5 


1.1 10~ 5 


1.3 


io- 


-2 


1.8 10~ 3 


7648 


2.4 


io- 


-3 


5.7 KT 4 


7.3 10" 


-4 


1.8 10~ 4 


1.2 10- 


-7 


1.5 10~ 8 


1.9 


io- 


-5 


2.7 10~ 6 


7560 


2.6 


10" 


-4 


3.0 10" 4 


2.7 10" 


-6 


3.1 10~ 6 


2.3 10" 


-4 


2.7 10" 4 


8.6 


10" 


-4 


9.9 10~ 4 


7549 


7.2 


io- 


-2 


4.1 10~ 2 


3.8 10" 


-2 


2.1 10~ 2 


1.6 10- 


-5 


8.1 10" 6 


1.8 


10" 


-3 


9.3 10~ 4 


7343 


3.0 


io- 


-2 


5.3 10- 3 


6.0 10" 


-2 


1.1 io- 2 


2.4 10- 


-5 


3.5 lO" 6 


4.2 


10" 


-3 


6.8 lO" 4 


6891 


2.3 


10" 


-2 


4.3 10~ 3 


1.2 10" 


-2 


2.2 10~ 3 


1.3 10" 


-6 


1.8 10" 7 


2.1 


10" 


-4 


3.2 10~ 5 


6707 


7.8 


10" 


-7 


7.4 10~ 7 


5.6 10" 


■21 


5.3 lO" 21 


1.7 10- 


-4 


1.5 10" 4 


6.7 


io- 


-4 


6.7 10" 4 


6533 


1.2 


io- 


-3 


3.7 10- 4 


1.6 10" 


-3 


5.1 10~ 4 


6.7 10- 


-7 


1.0 10- 7 


1.2 


io- 


-4 


2.1 10" 5 


6225 


6.7 


io- 


-4 


2.8 10" 4 


1.5 10" 


-4 


6.2 lO" 5 


1.0 io- 


-8 


1.8 lO" 9 


2.0 


10" 


-6 


3.9 lO" 7 
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Table 3: Characteristics for a subsample of long duration GRBs extracted from the sample 
of 566 long duration GRBs. The luminoisty of each burst is calculated from its variability, as 
described in the text. The fluence and zenith angle of each burst are taken from the BATSE 
catalog, wheras the redshift is derived from the luminosity and fluence. The break energy 
is obtained by fitting the BATSE data, and the Lorentz factor is calculated as described in 
Appendix A. 



BATSE # 


F 7 ( 


'erg 


;/cm 2 ) 


L 7 (erg/s) 




z 


£5(MeV) 


r 


9 Z (degrees) 


676 


4. 


19 


10" 


-5 


1.7 


10 52 


1 


.883 


0.260 


273. 


135.2 


1606 


6 


,20 


10" 


-5 


1.3 


10 52 


1 


.313 


0.280 


261. 


45.2 


2102 


2. 


,47 


10" 


-6 


4.6 


1Q 50 





.644 


0.330 


150. 


34.7 


2431 


1 


,84 


io- 


-5 


1.1 


10 51 





.251 


0.480 


173. 


71.0 


2586 


1 


.45 


io- 


-5 


2.9 


1Q 53 


4 


.487 


0.500 


437. 


99.3 


2798 


2. 


.31 


io- 


-4 


1.9 


10 51 





.381 


0.460 


190. 


30.0 


3356 


1 


.60 


10" 


-6 


7.9 


10 51 


1 


.837 


0.230 


241. 


65.8 


5644 


1 


.22 


io- 


-5 


9.5 


10 52 


3 


.376 


0.380 


364. 


131.9 


6397 


2 


.11 


10" 


-5 


9.3 


10 51 


1 


.306 


0.490 


247. 


81.6 


6672 


7. 


.81 


io- 


-6 


1.6 


10 53 


4 


.311 


0.250 


398. 


108.3 


7822 


9 


,09 


io- 


-6 


2.3 


10 52 


2 


.417 


0.420 


287. 


134.0 


8008 


1 


.07 


10- 


-4 


2.8 


10 50 





.255 


1.940 


138. 


144.3 
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Table 4: Estimated neutrino event rates for the sample of long duration GRBs described in 

table 3. Rates for muon and shower events are shown for the IceCube neutrino telescope. 

Tau events are not shown. Note that in the supranova scenario (models 3 or 4), internal 

shock processes (models 1 or 2) also contribute to the event rate. The total predicted rate 

for models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the rates 

shown for the internal shock model. 

Model 1 Model 2 Model 3 Model 4 



BATSE # 








Shower 








Shower 




H 




Shower 




n 




Shower 


676 


3.1 


10" 


-3 


6.9 


10" 4 


8.5 


io- 


-3 


1.9 IO" 3 


3.6 


10" 


-6 


5.1 


io- 7 


6.1 


10" 


-4 


9.7 10^ 5 


1606 


4.8 


10" 


-3 


4.1 


io- 3 


1.4 


io- 


-2 


1.2 IO" 2 


6.6 


io- 


-6 


5.4 


10~ 6 


8.3 


io- 


-4 


6.8 HT 4 


2102 


2.6 


io- 


-4 


2.6 


10" 4 


3.8 


10" 


-4 


3.7 10" 4 


5.5 


10" 


-8 


5.3 


io- 8 


7.3 


io- 


-6 


7.0 10~ 6 


2431 


1.6 


10" 


-3 


6.0 


10~ 4 


2.6 


io- 


-3 


1.0 10~ 3 


1.3 


10" 


-6 


4.9 


io- 7 


1.5 


10" 


-4 


5.8 10~ 5 


2586 


2.9 


io- 


-3 


6.5 


io- 4 


6.3 


10" 


-3 


1.4 10~ 3 


3.1 


10" 


-5 


4.1 


10~ 6 


3.1 


10" 


-2 


3.9 IO" 3 


2798 


1.9 


io- 


-2 


2.0 


io- 2 


3.5 


10" 


-2 


3.7 IO" 2 


8.1 


10" 


-4 


8.2 


io- 6 


9.4 


io- 


-4 


9.4 IO" 4 


3405 


1.3 


10" 


-3 


1.1 


IO" 3 


2.5 


io- 


-3 


2.2 10~ 3 


1.7 


10" 


-5 


1.5 


IO" 5 


1.3 


10" 


-2 


1.1 IO" 2 


5644 


1.4 


io- 


-3 


3.5 


10~ 4 


3.3 


10" 


-3 


8.3 IO" 4 


3.1 


io- 


-6 


4.6 


io- 7 


4.8 


io- 


-4 


8.2 IO" 5 


6397 


4.3 


io- 


-3 


8.1 


10~ 4 


7.6 


10" 


-3 


1.5 IO" 3 


9.0 


io- 


-6 


1.6 


IO" 6 


1.1 


10" 


-3 


2.0 10~ 4 


6672 


1.9 


io- 


-3 


2.8 


10" 4 


6.0 


10" 


-3 


8.6 10" 4 


8.4 


10" 


-6 


1.1 


io- 6 


9.9 


io- 


-4 


1.6 KT 4 


7822 


9.7 


10" 


-4 


1.7 


10~ 4 


1.8 


io- 


-3 


3.2 10~ 4 


1.0 


10" 


-6 


1.5 


io- 7 


1.8 


10" 


-4 


2.9 IO" 5 


8008 


2.4 


10" 


-2 


8.7 


10~ 3 


7.4 


io- 


-3 


2.7 10~ 3 


8.4 


10" 


-7 


1.1 


io- 7 


1.7 


io- 


-4 


2.5 IO" 5 
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Table 5: Characteristics for a subsample of short duration GRBs extracted from the larger 
sample of 199 short duration GRBs. The redshift of each event is fixed to z — 1. The fluence 
and zenith angle of each burst are taken from the BATSE catalog, whereas the luminosity is 
derived from the redshift and fluence. The break energy is obtained by fitting the BATSE 
data, and the Lorentz factor is calculated as described in the Appendix A. 



BATSE # 


F 7 [ 


'erg 


l/cm 2 ) 


L 7 (erg/s) 




£?(MeV) 


r 


9 Z (degrees) 


603 


8. 


.93 


10" 


-7 


1.5 


10 52 


1 


0.480 


392. 


54.2 


1073 


2. 


,89 


10" 


-7 


2.4 


10 52 


1 


0.490 


425. 


101.2 


1518 


1 


.32 


10" 


-6 


7.4 


10 51 


1 


0.240 


1558. 


41.3 


2312 


5. 


.75 


10" 


-7 


2.6 


10 52 


1 


0.280 


431. 


72.2 


2861 


1 


.77 


10" 


-6 


4.0 


10 51 


1 


0.390 


315. 


130.8 


3215 


1 


.38 


10" 


-6 


4.5 


10 52 


1 


0.270 


471. 


106.6 


3940 


3 


.66 


io- 


-7 


1.2 


10 52 


1 


0.080 


2390. 


71.1 


5212 


7. 


.35 


10" 


-7 


2.5 


10 52 


1 


0.270 


427. 


28.0 


6299 


2. 


,14 


10" 


-7 


3.6 


10 52 


1 


0.180 


1208. 


35.7 


7063 


1 


,41 


10" 


-6 


1.0 


1Q 53 


1 


0.280 


541. 


36.3 


7447 


2 


.86 


10" 


-7 


2.5 


10 52 


1 


0.790 


427. 


63.1 


7980 


2 


,15 


10" 


-7 


3.2 


1Q 52 


1 


0.070 


2006. 


1.6 
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Table 6: Estimated neutrino event rates for the sample of short duration GRBs described in 

table 5. Rates for muon and shower events are shown for the IceCube neutrino telescope. 

Tau events are not shown. Note that in the supranova scenario (models 3 or 4), internal 

shock processes (models 1 or 2) also contribute to the event rate. The total predicted rate 

for models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the rates 

shown for the internal shock model. 

Model 1 Model 2 Model 3 Model 4 



BATSE # 




fi 




Shower 








Shower 










Shower 




fi 




Shower 


603 


6.2 


10' 


-5 


4.4 


IO" 5 


2.3 


io- 


-4 


1.6 10" 4 


1 


.9 


io- 


-8 


1.3 10~ 8 


2.6 


io- 


-6 


1.8 IO" 6 


1073 


2.3 


10" 


-5 


3.8 


10~ 6 


9.0 


io- 


-5 


1.5 10~ 5 


2 


.0 


io- 


-8 


2.6 10~ 9 


3.0 


io- 


-6 


4.4 10~ 7 


1518 


5.6 


io- 


-6 


4.9 


IO" 6 


1.5 


10" 


-7 


1.3 10" 7 


1 


.5 


io- 


-6 


1.3 10~ 6 


5.1 


io- 


-5 


4.5 KT 5 


2312 


2.3 


io- 


-5 


8.4 


10~ 6 


1.1 


io- 


-4 


3.9 10" 5 


3 


.0 


io- 


-8 


1.1 10~ 8 


4.0 


io- 


-6 


1.4 KT 6 


2861 


5.1 


io- 


-5 


8.3 


10~ 6 


1.7 


io- 


-4 


2.7 10~ 5 


1. 


.5 


io- 


-8 


2.2 10~ 9 


2.9 


io- 


-6 


4.8 KT 7 


3215 


7.5 


io- 


-5 


1.3 


10~ 5 


3.6 


io- 


-4 


6.2 IO" 5 


9 


.1 


io- 


-8 


1.2 IO" 8 


1.4 


io- 


-5 


2.1 IO" 6 


3940 


2.9 


io- 


-6 


1.1 


10~ 6 


6.6 


io- 


-8 


2.6 10~ 8 


6 


.8 


io- 


-6 


2.6 IO" 6 


3.6 


io- 


-5 


1.4 KT 5 


5212 


2.3 


io- 


-5 


2.4 


10~ 5 


1.1 


io- 


-4 


1.1 10' 4 


1. 


.3 


io- 


-8 


1.3 IO" 8 


1.8 


io- 


-6 


1.8 IO" 6 


6299 


1.6 


io- 


-6 


1.6 


10~ 6 


7.4 


io- 


-7 


7.0 IO" 7 


1 


.0 


io- 


-7 


9.9 10" 8 


1.3 


io- 


-5 


1.2 10~ 5 


7063 


8.5 


io- 


-6 


8.0 


10~ 6 


4.2 


10" 


-5 


4.0 10" 5 


5. 


.0 


io- 


-8 


4.7 10~ 8 


5.8 


io- 


-6 


5.5 10~ 6 


7447 


3.7 


io- 


-5 


2.2 


10~ 5 


1.1 


10" 


-4 


6.6 10" 5 


9 


.9 


io- 


-9 


5.2 IO" 9 


1.3 


io- 


-6 


7.0 IO" 7 


7980 


3.6 


io- 


-6 


4.4 


io- 6 


5.0 


io- 


-7 


6.0 10' 7 


3 


.1 


io- 


-7 


3.7 IO" 7 


6.2 


io- 


-6 


7.2 10~ 6 
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Table 7: Characteristics for X-ray flash candidates. The redshift to each event is fixed to 
z — 1. The fluence and zenith angle of each burst are taken from the BATSE catalog, whereas 
the luminosity is derived from the redshift and fluence. The break energy is obtained by 
fitting the BATSE data, and the Lorentz factor is calculated as described in Appendix A. 



BATSE # 


F 7 (erj 


;/cm 2 ) 


L 7 (erg/s) 


z 


£ 7 (MeV) 


r 


9 Z (degrees) 


659 


7.60 


10~ 6 


9.9 10 48 


1 


0.140 


1508. 


109.4 


717 


5.44 


io- 7 


1.0 10 49 


1 


0.050 


2505. 


98.1 


927 


6.31 


io- 7 


1.7 10 49 


1 


0.130 


1357. 


128.4 


1244 


3.12 


io- 6 


1.9 10 49 


1 


0.060 


1961. 


133.7 


2381 


4.96 


io- 7 


1.1 10 50 


1 


0.050 


1367. 


79.8 


2917 


5.27 


10~ 7 


3.8 10 49 


1 


0.160 


1004. 


94.6 


2990 


4.78 


io- 7 


5.3 10 49 


1 


0.050 


1654. 


101.6 


3141 


1.55 


io- 7 


5.4 10 50 


1 


0.050 


926. 


141.7 


5483 


4.18 


io- 7 


6.4 10 50 


1 


0.090 


663. 


154.8 


5497 


3.21 


10~ 6 


7.1 10 48 


1 


0.060 


2500. 


141.1 


5627 


3.16 


io- 7 


8.4 10 49 


1 


0.070 


1247. 


63.4 


6167 


9.21 


10~ 6 


2.3 10 49 


1 


0.070 


1724. 


126.3 


6437 


4.56 


io- 7 


1.2 10 50 


1 


0.050 


1359. 


118.4 


6582 


1.17 


io- 6 


4.8 10 50 


1 


0.060 


872. 


27.7 


7942 


9.86 


io- 7 


4.5 10 49 


1 


0.050 


1722. 


134.5 
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Table 8: Estimated neutrino event rates for X-ray flash candidates. Rates for muon and 

shower events are shown for the IceCube neutrino telescope. Tau events are not shown. 

Note that in the supranova scenario (models 3 or 4), internal shock processes (models 1 or 

2) also contribute to the event rate. The total predicted rate for models 3 or 4 is, therefore, 

the sum of rates shown for the supranova model and the rates shown for the internal shock 

model. 

Model 1 Model 2 Model 3 Model 4 



BATSE # 








Shower 


n 




Shower 








Shower 








Shower 


659 


4.8 


io- 


-5 


4.6 


10" 


-6 


6.6 10- 


ii 


6.3 10" 


-12 


1.5 


io- 


-4 


1.7 10" 


-5 


2.3 


io- 


-3 


5.1 10" 


-4 


717 


1.3 


10" 


-6 


1.7 


io- 


-8 


6.6 10- 


-12 


8.6 10- 


-15 


4.1 


io- 


-5 


4.8 10- 


-6 


2.3 


io- 


-4 


5.7 10- 


-5 


927 


2.3 


io- 


-5 


6.8 


io- 


-6 


9.1 10" 


-11 


2.7 10- 


-11 


3.6 


io- 


-6 


3.9 10- 


-7 


1.2 


io- 


-4 


2.8 10- 


-5 


1244 


9.1 


io- 


-7 


7.0 


io- 


-8 


1.9 10- 


-13 


1.5 10- 


-13 


2.5 


io- 


-5 


3.1 10- 


-6 


6.6 


io- 


-4 


1.8 10- 


-4 


2381 


1.1 


io- 


-5 


2.4 


io- 


-6 


7.3 10" 


-10 


1.6 10" 


-10 


2.4 


io- 


-5 


5.1 10- 


-6 


1.8 


io- 


-4 


4.4 10" 


-5 


2917 


8.6 


io- 


-5 


1.8 


io- 


-5 


2.0 10 


-9 


4.3 10- 


-10 


2.4 


io- 


-5 


2.3 10- 


-6 


1.9 


io- 


-4 


3.61 0- 


-5 


2990 


1.4 


io- 


-6 


1.1 


io- 


-7 


2.0 10- 


-11 


1.5 10- 


-12 


1.8 


io- 


-5 


2.0 10- 


-6 


1.7 


io- 


-4 


3.8 10- 


-5 


3141 


4.0 


io- 


-7 


7.3 


io- 


-8 


5.9 10- 


-10 


1.1 io- 


-10 


2.6 


io- 


-7 


2.7 10- 


-8 


1.7 


io- 


-5 


3.7 10- 


-6 


5483 


6.7 


io- 


-7 


7.4 


io- 


-8 


2.5 10 


-9 


2.7 10- 


-10 


2.5 


io- 


-7 


2.3 10- 


-8 


1.8 


io- 


-5 


3.9 10- 


-6 


5497 


4.5 


io- 


-7 


3.5 


io- 


-8 


1.4 10" 


-13 


1.0 10" 


-14 


2.8 


io- 


-5 


3.9 10- 


-6 


6.8 


io- 


-4 


2.1 10- 


-4 


5627 


4.9 


io- 


-6 


2.6 


io- 


-6 


2.4 10- 


-10 


1.3 10- 


-10 


5.8 


io- 


-6 


3.0 10" 


-6 


5.1 


io- 


-5 


2.7 10- 


-5 


6167 


5.4 


io- 


-6 


4.1 


io- 


-7 


2.0 10" 


-11 


1.5 10" 


-12 


8.1 


10" 


-5 


9.6 10- 


-6 


2.0 


io- 


-3 


5.2 10" 


-4 


6437 


4.5 


io- 


-7 


3.5 


10" 


-8 


3.1 10" 


-11 


2.4 10- 


-12 


4.2 


10" 


-6 


4.6 10- 


-7 


1.0 


io- 


-4 


2.4 10" 


-5 


6582 


1.2 


io- 


-5 


1.2 


io- 


-5 


1.6 10 


-8 


1.7 10 


-8 


8.5 


io- 


-6 


8.9 10- 


-6 


8.0 


io- 


-5 


8.2 10" 


-5 


7942 


2.6 


io- 


-7 


1.7 


io- 


-8 


2.7 lO- 


-12 


1.7 10- 


-13 


6.3 


io- 


-6 


7.5 10- 


-7 


1.9 


io- 


-4 


5.0 10- 


-5 



Table 9: Sum of event rates and event rate per GRB in a kilometer scale neutrino telescope 
for GRBs with known redshift. Note that in the supranova scenario (models 3 or 4), internal 
shock processes (models 1 or 2) also contribute to the event rate. The total predicted rate 
for models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the rates 
shown for the internal shock model. 



Model 


total fi 


total showers 


total r 


/i/GRB 


showers/ GRB 


r/GRB 


1 


0.236 


7.49 10" 2 


2.05 10" 3 


1.57 10~ 2 


4.99 10" 3 


1.37 10" 4 


2 


0.215 


5.72 10~ 2 


1.26 10~ 3 


1.44 10~ 2 


3.81 10~ 3 


8.43 10~ 5 


3 


5.30 10" 4 


4.45 10~ 4 


2.50 10- 4 


3.54 10" 5 


2.97 10~ 5 


1.67 10- 5 


4 


2.12 10~ 2 


5.19 10- 3 


4.57 10" 4 


1.41 10~ 3 


3.46 10~ 4 


3.05 10~ 5 
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Table 10: Sum of event rates and event rate per GRB in a kilometer scale neutrino telescope 
for long duration GRBs. Note that in the supranova scenario (models 3 or 4), internal shock 
processes (models 1 or 2) also contribute to the event rate. The total predicted rate for 
models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the rates 
shown for the internal shock model. 



Model 


total /i 


total showers 


total r 




showers/ GRB 


r/GRB 


1 


3.54 


1.09 


5.18 10" 2 


6.16 10" 3 


1.90 10" 3 


9.02 10" 5 


2 


5.07 


1.65 


4.76 10~ 2 


8.82 10- 3 


2.87 10~ 3 


8.29 10~ 5 


3 


4.26 10~ 2 


1.22 10~ 2 


1.41 10~ 2 


7.41 10~ 5 


2.12 10~ 5 


2.45 10~ 5 


4 


3.18 


0.87 


1.27 10- 1 


5.53 10~ 3 


1.51 10~ 3 


2.20 10" 4 



Table 11: Sum of estimated event rates and event rate per GRB in a kilometer scale neutrino 
telescope for short duration GRBs. Note that in the supranova scenario (models 3 or 4), 
internal shock processes (models 1 or 2) also contribute to the event rate. The total predicted 
rate for models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the 
rates shown for the internal shock model. 



Model 


total \x 


total showers 


total r 


/i/ GRB 


showers/ GRB 


r/GRB 


1 


2.17 10" 2 


6.73 10" 3 


2.41 10" 4 


1.09 10" 4 


3.38 10" 5 


1.21 10" 6 


2 


8.39 10~ 2 


2.58 10~ 2 


8.85 10~ 4 


4.22 10~ 4 


1.30 10~ 4 


4.45 10~ 6 


3 


1.18 10~ 4 


4.81 10~ 5 


4.42 lO" 5 


5.94 10~ 7 


2.42 10~ 7 


2.22 10~ 7 


4 


7.69 10~ 3 


2.04 10- 3 


6.60 10" 4 


3.86 10~ 5 


1.02 10~ 5 


3.32 10~ 6 



Table 12: Sum of estimated event rates and event rate per GRB in a kilometer scale neutrino 
telescopes for X-ray flash candidates. Note that in the supranova scenario (models 3 or 4), 
internal shock processes (models 1 or 2) also contribute to the event rate. The total predicted 
rate for models 3 or 4 is, therefore, the sum of rates shown for the supranova model and the 
rates shown for the internal shock model. 



Model 


total /i 


total showers 


total r 


yu/GRB 


showers/ GRB 


r/GRB 


1 


1.96 10- 4 


4.74 10~ 5 


3.24 10" 5 


1.30 10" 5 


3.16 10" 6 


2.16 10" 6 


2 


2.22 lO" 8 


1.79 10~ 8 


1.59 10~ 8 


1.48 10~ 9 


1.19 10~ 9 


1.06 10" 9 


3 


4.23 10~ 4 


6.09 10~ 5 


2.51 10~ 5 


2.82 10~ 5 


4.06 10~ 6 


1.67 10~ 6 


4 


7.00 10~ 3 


1.80 10~ 3 


2.14 10~ 5 


4.67 10~ 4 


1.20 10~ 4 


1.43 10~ 6 
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Fig. 1. — The distribution of (a) redshifts, (b) Lorentz factors, (c) ffuences and (d) zenith 
angles (for muon events: dotted line corresponds to model 1, solid line to model 2) for long 
duration GRBs. In (a), the decrease at large redshift is due to sampling bias. The Lorentz 
factors in (b) were calculated as described by Eqs. A17 and A18. Frame (d) demonstrates 
the advantages of long muon range and poor absorption near the horizon for muon track 
detection. 
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Fig. 2. — The distribution of the estimated number of muon events from individual (a) long 
duration GRBs (the dotted line corresponds to model 1 and the solid line to model 2) and 
(b) X-ray flash candidates (dotted line corresponds to model 3 and the solid line to model 
4). Note that the majority of events result from a relatively small number of GRBs. 
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Fig. 3. — The muon neutrino spectrum, e v {dN v l de v \ for our fiducial parameters in models 1 
and 2 (interactions with GRB photons). The solid line is for a typical GRB with T = 300 and 
t v = 10 ms, while the dashed line is for a X-ray flash candidate with T = 1440 (calculated 
from Eq. A18) and t v = 50 ms. 
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Fig. 4. — The proton energies that correspond to break energies in the neutrino spectrum, 
£k ~ as a function of t S( j. This figure is taken from Guetta & Granot 2002a 
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Fig. 5. — The muon neutrino spectrum, el(dN u /de u ), for our fiducial parameters in models 
3 and 4 (interactions with PWB photons). Four choices of t s d = 0.01, 0.07, 0.4, 30 yr are 
used, which correspond to the 4 different orderings of the break energies. This figure is taken 
from Guetta & Granot 2002a 
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Fig. 6. — The probabilities of a neutrino generating various types of events when traveling 
through the effective area of a neutrino telescope. Curves are shown for various choices of 
zenith angle, which reflects enhancements due to long muon range and the effect of absorption 
in the Earth. 
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Fig. 7. — The probability of a muon neutrino generating a detectable muon as a function of 
zenith angle. The curves show the probabilities for three representative neutrino energies: 
10 12 eV, 10 15 eV and 10 18 eV. 



